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Abstract 

P 

We evaluate the threshold matrix-element for the reaction pp — > pp-yp in a fully relativis- 
tic Feynman diagrammatic approach. We employ a simple effective range approximation 
to take care of the S-wave pp final-state interaction. The experimental value for the 
threshold amplitude A = (2.7 — i 0.3) fm 4 can be reproduced by contributions from 
tree level chiral (long-range) pion exchange and short-range effects related to vector 
meson exchanges, with u exchange giving the largest individual contribution. Pion loop 
effects appear to be small. We stress that the commonly used heavy baryon formalism 
is not applicable in the NN-system above the pion production threshold due to the 
large external momentum, \p\ ~ yJMm^, with M and m n the nucleon and the pion 
mass, respectively. We furthermore investigate the reaction pp —> pnn + near threshold 
within the same approach. We extract from the data the triplet threshold amplitude 
as B = (2.8 — il.5) fm 4 . Its real part can be well understood from (relativistic) tree 
level meson-exchange diagrams. In addition, we investigate the process pp — > ppn near 
threshold. We use a simple factorization ansatz for the ppn final-state interaction and 
extract from the data the modulus of the threshold amplitude, \C\ = 1.32 fm 4 . With 
9r)N = 5.3, this value can be reproduced by (relativistic) tree level meson-exchange dia- 
grams and 77-rescattering, whose strength is fixed by the nN scattering length. We also 
comment on the recent near threshold data for //-production. 
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1 Introduction and summary 



With the advent of proton cooler synchrotrons at Bloomington, Jiilich and Uppsala, high 
precision data for the processes pp — > ppir , pp — > dir + , pp — > pmr + and pp — > pprj in the 
threshold region [I], [^, [|, have become available. The first data on neutral pion production 
were a big surprise because the experimental cross sections turned out to be a factor of five 
larger than the theoretical predictions based on direct pion production and neutral pion 
rescattering fixed from on-shell ttN data. Nevertheless, the energy-dependence of the total 
cross section cr tot from threshold to excess energies of about 50 MeV is completely given by 
the strong S-wave pp final-state interaction, see e.g. ref.@. Subsequently, it was argued that 
heavy-meson exchanges might be able to remove this discrepancy [0, E|. On the other hand 
it was shown |10j that the (model-dependent) off-shell behavior of the full nN T-matrix 
also can enhance the cross section considerably. Another avenue to eventually understand 
the near threshold cross section is offered by calculations within the framework of tree-level 



heavy baryon chiral perturbation theory (HBChPT) including dimension two operators |IT 
|l2t p~3| , |i~4jl . In these papers, HBChPT has been used to constrain the long-range 7r°-exchange 
contributions. However, it was found that the calculations for pp — > ppn performed so far 
lead to a marked difference in the role of the so-called 7r° rescattering contribution, which 
interferes constructively with the direct production in the Jiilich model and destructively in 
the HBChPT framework, respectively. In a recent paper, it was argued that the treatment 
underlying the isoscalar pion-nucleon scattering amplitude and the related transition operator 
for the process NN —>■ NNn (where N denotes the nucleon) in the HBChPT framework is not 
yet sufficiently accurate and thus the resulting rescattering contribution should be considered 
as an artifact of this approximation [T3j . It was also stressed that the process pp — > dix + is more 
sensitive to the long-range (chiral) pion exchange. One-loop graphs have been considered 
in 
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However, one important feature of this reaction should be kept in mind, namely the large 
momentum mismatch between the initial and the final nucleon-nucleon state: the initial 
center-of-mass momentum is given by \p | 2 = Mm^ + with (M) the neutral pion 

(proton) mass, whereas close to threshold the final one is compatible with zero. This large 
difference leads to momentum transfers equal to \p\ and higher. In particular, both the 
squared invariant momentum transfer and the kinetic energy of the incoming protons are 
—Mm n and respectively, and therefore the chiral counting in the original sense does 

not apply, see e.g. (T^|. We will actually demonstrate here that non-relativistic approximations 



(i.e. the heavy baryon formalism or HBChPT) are not applicable in the two-nucleon system 
above the pion production threshold. The simple reason for that is the extreme kinematics 
with the external momentum \p\ ~ \jMm^ diverging in the heavy nucleon limit M — > oo 
(while this argument could be considered formal, it nicely pinpoints the problems of the heavy 
baryon approach). The HBChPT framework (which is very successful for low-energy single 
nucleon processes) looses its systematic order-by-order low-energy expansion when external 
momenta grow with the nucleon mass (\p\ ~ \M). In order to avoid such problems (which are 
simply related to kinematics) it seems mandatory to perform fully relativistic calculations. In 
particular, we find that one of the low-energy constants C2 entering the chiral 7r°-rescattering 
is enhanced in the relativistic calculation by a factor of two as compared to the non-relativistic 
HBChPT approach. In this spirit we perform here a fully relativistic calculation of various 
tree and one-loop diagrams contributing to pp — > pp7T° and pp —>■ pmr + at threshold. Our 
approach can also be extended in a straightforward manner to r] and rj' production in pp- 
collisions. We thus discuss here also the reaction pp — > pprj up to laboratory excess energies 
of about 100 MeV using a simple factorized form to account for the final-state interactions 
in the pprj three-body system. 
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The pertinent results of this investigation can be summarized as follows: 

(i) Assuming that the pp final-state interaction is an on-shell NN-process and using a 
simple effective range parametrization for the 1 5'o pp phase shift, we can accurately fit 
the 40 data points of the total cross section from threshold up to Ti ab = 326 MeV with 
a constant (threshold) amplitude equal to A — (2.7 — i 0.3) fm 4 . 

(ii) The real part of this number can be well understood in terms of chiral ir° exchange 
(including chiral tt° rescattering) and heavy meson (u, p°, rj) exchanges based on a 
relativistic Feynman diagram calculation. 

(iii) We have evaluated some classes of one-loop graphs and find that they lead to small 
corrections of the order of a few percent. Therefore chiral loops do not seem to play 
any significant role in the processes NN — > NNn, which are dominated by one-pion 
exchange and short-range physics (with the notable exception of pp — > dn + , where d 
denotes the deuteron). 

(iv) Both the long range 7r° exchange and the short range vector meson exchange lead to 
contributions to the threshold amplitude A which do not vanish in the chiral limit 
m n — > 0. There is no chiral suppression of the reaction pp — > ppir compared to other 
NNn channels. In all cases the respective threshold amplitudes are non-zero (and 
finite) in the chiral limit. This is in contrast to the widespread believe that pp — > ppn° 
is suppressed for reasons of chiral symmetry. 

(v) Within the same approach, we have investigated the threshold behavior of the process 
pp — > pri7[ + . It is given in terms of A and the triplet threshold amplitude B with 
the empirical value B = (2.8 — % 1.5) fm 4 . The corresponding real part Re£> is well 
reproduced by chiral one-pion exchange and short-range vector meson physics. The 
empirical value of B has however a sizeable imaginary part, which naturally can not 
be explained by tree graphs. This channel deserves some further study. The small and 
large imaginary parts of A and B, respectively, are related to the weak and strong initial 
state interaction in the 3 Po an d 3 -Pi entrance channel, respectively. 

(vi) Finally, we have investigated the process pp — > pprj near threshold. We treat the final 
state-interaction in the pp and rjp subsystems using a factorization ansatz and effective 
range approximations. The empirical value of the modulus of the threshold amplitude, 
\C\ = 1.32 fm 4 , is reproduced with = 5.3 by tree level meson exchange graphs 
and 77-rescattering. The strength of the latter is fixed from the (real part of the) i]N 
scattering length. 

(vii) Altogether, it is quite surprising that the total cross sections for pion and eta produc- 
tion in pp-collisions near threshold can be explained so simply in terms of final-state 
interaction (using effective range approximations) and certain well-known (relativistic) 
meson-exchange diagrams for the respective threshold amplitudes. However, such an 
approach requires further improvements and tests. First, one should include system- 
atically initial-state and final-state NN-interactions. Secondly, in order to perform a 
detailed test of the underlying meson production mechanism one has to consider more 
exclusive observables like angular distributions of differential cross sections and asym- 
metries generated by polarized proton beams and targets. We hope to report on these 
topics in the near future. 
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2 Threshold kinematics and final— state interaction 



2.1 Kinematics 

We consider the reaction pi(p) + P2{—p) — > P + V + tt° in the center-of-mass (cm) frame at 
threshold (see fig. [I]). The invariant T-matrix can be expressed in terms of one complex- 
valued (constant) amplitude, which we denote by A, as 

T^(pp^ppn°)=A(ia 1 -ia 2 + a 1 xa 2 )-p . (1) 

The <7i,2 are the spin-matrices of the two protons. The value of the proton cm momentum to 
produce a neutral pion at rest is given by 

\p | = y/m„(M + m n /4) = 362.2 MeV , (2) 

with M = 938.27 MeV the proton and m n = 134.97 MeV the neutral pion mass, respectively. 
Obviously, \p\ vanishes in the chiral limit of zero pion mass. Therefore the soft-pion theorem 
which requires a vanishing threshold T-matrix in the chiral limit = is trivially fulfilled 
(as long as A does not become singular) [|17]]. We remark that similar features occur for the 



reaction nN — > nnN (see ref.[|T8|). All dynamical information is encoded in the threshold 
amplitude A of dimension [length 4 ]. In the threshold region, the wave function of the final 
di-proton system as well as that of the neutral pion are dominated by angular momentum 
zero states, thus we are dealing with a 3 P — > 1 5 , o s transition. Consequently, one deduces 
from unitarity 

A = \A\ e lS{3p ^ , (3) 

with the 3 Po PP phase shift to be taken at the threshold energy in the lab frame, T^ 3 = 
m^(2 + m n /2M) = 279.65 MeV, where 5{ 3 P ) = -6.3° (FA95 solution of VPI). Thus the 
imaginary part Im^4 is about —1/9 of the real part Re A and contributes negligibly to the total 
cross section near threshold proportional to \A\ 2 . The threshold T-matrix is a pseudoscalar, 
it is symmetric under the exchange of the two ingoing protons d\ <-> a 2 , p — > —p- Furthermore 
Eq. ([!]) incorporates the Pauli exclusion principle for the (indistinguishable) outgoing protons, 
since left multiplication with the spin exchange operator (l + c?i • <?2)/2 leads to a minus sign 
by the identity: |(1 + B\ ■ cr 2 ) (i&i — i &2, &i x &i) = — (&i x a 2 , i (Ti — i a 2 )- Diagrams with 
crossed proton lines are therefore automatically included. Approximating the near threshold 
T-matrix by the T-matrix exactly at threshold one gets for the unpolarized total cross section 



^tot(Tlab) - 167r 2^ 2 + ^)9/2 l^' 2 ( Tlab _ T ^) 2 > ( 4 ) 

with /i = m-x/M. Note that the flux and three-body phase space factors have been approx- 
imated in Eq.(|J) by an analytical expression which is accurate within a few percent in the 
threshold region. Such a form has already been proven to be quite accurate in chiral pertur- 
bation theory studies of the reactions 7p — > tt tc°p \ 1J| and tcN — > tcttN [ 20| , where 7 denotes 
a real photon. In the case of pp — > ppn , however, such an approximation is not sufficient near 
threshold. This can be seen e.g. by taking the near threshold data of refs. |], [J and dividing 
them by the three-body phase space factor ~ (T lab - T^) 2 , i.e. a to t(Piab) = C ■ (T lab - T^) 2 . 
The resulting values for C are not constant in energy. This, of course, has to do with the 
strong pp final-state interaction in the 1 So partial wave. So before we can extract a value for 
we have to correct for the final-state interaction. 
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2.2 Treatment of the final— state interaction 

The final-state interaction in the 1 S'o di-proton state modifies the simple phase-space formula 



Eq.([|). We follow here a procedure derived by Watson [51] where one essentially assumes that 
final-state interaction is taking place only when the nucleons are on their mass-shell. In this 
approach, the unpolarized total cross section for pp — > ppir including final-state interaction 
takes the form 



(T \ I A \2 f M \ 3 2yXiab 

^totUlabj - \A\ { — ) (2M + Tiab)3/2 



/•Wmax / 

x dWJ(W 2 -4M 2 )\(W 2 ,ml,4M 2 + 2MT lah )F p (W) . (5) 

Jim v 

Here, F p (W) is the correction factor due the final-state interaction. We evaluate it in the 
effective range approximation. This is of course a very strong assumption but it allows to 
explain the energy dependence of the experimental total cross sections very accurately in 
terms of a single constant amplitude A. Separating off the final-state interaction in that way, 
we can then pursue a diagrammatic approach to the (on-shell) production amplitude A. This 
will allow us to investigate in a simple fashion the role of one-pion exchange and chiral loop 
effects together with shorter range exchanges due to heavier mesons. For a more complete 
dynamical description including off-shell final-state interactions and so on, a microscopic 
model treating also (unobservable) half-off shell effects is needed. After these remarks, we 
return to the effective range approximation, in which the function F P (W) takes the form 

W = ff'-S) = I 1 + f <°> + ^ - 4M2 » + % w2 ' 4M2 > 2 P < 6 > 

with W the final-state di-proton invariant mass and X(x,y,z) = x 2 + y 2 + z 2 — 2yz — 2xz — 2xy 
the Kallen function. W max = \JAM 2 + 2MTi a b — is the kinematical endpoint of the 
di-proton invariant mass spectrum. Furthermore, a p = (7.8098 ± 0.0023) fm and r p — 



[2.7Q7 ± 0.010) fm, taken from ref.[^], are the scattering length and effective range pa- 
rameter for elastic pp-scattering including electromagnetic effects. Note that we have fixed 
the normalization of the correction factor F P (W) such that in the limit of vanishing scatter- 
ing length a p (i.e. vanishing final-state interaction), F p (W) becomes identical to one (in the 
effective range approximation). Furthermore, the condition F p (2M) = 1 ensures that there is 
no final-state interaction effect exactly at threshold, as it must be according to the definition 
of the threshold amplitude A in Eq.(l). In appendix C we give a simple derivation of F p (W) 
in scattering length approximation (i.e. for r p = 0) using effective field theory methods. 

With the help of Eqs.(|,^]) we are now in the position to extract the values of |^4| as shown 
in table We ignore the data point at the lowest energy since that close to threshold a better 
treatment of the infinite range Coulomb interaction in terms of proton wave functions would 
be needed. The empirical value of A is thus 

A {cxp) = (2.7 -i 0.3) fm 4 , (7) 

anticipating the positive sign of the real part from the calculation in section 3. Since we are 
mostly interested in achieving a qualitative picture of the underlying production process, we 
refrain from assigning an uncertainty to this number. This number should be comparable to 
the one found by Adam et al. [p3fl , who employ essentially the same method to correct for 



the final-state interaction but use a different (unspecified) normalization of the T-matrix. In 



Fig. 



we show the resulting total cross sections in comparison to the data of refs. |], using 



Eqs.(|J|,0). 
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Ti ab [MeV] 


282.2 285 290 295 300.3 308 314 319 325.6 


ftot [/-tb] 


0.148 0.56 1.31 2.06 3.07 4.47 5.25 6.18 7.71 


|^| [fm 4 ] 


2.42 2.73 2.69 2.64 2.70 2.71 2.64 2.66 2.73 



Table 1: Extracted values for the threshold amplitude |^4| for the proton laboratory kinetic 
energies Ti a b and cross sections er t ot of ref.0. 



3 Meson— exchange contributions 

In this section, we work out a variety of one-boson exchange contributions related to chiral 
(pion) and non-chiral (vector meson) physics to the amplitude A. 

3.1 Tree level Goldstone boson contributions 

In this and the following paragraph, we consider tree level and one-loop pion exchanges 
contributing to the ppn production process at threshold. Due to the large proton momentum 
even at threshold, \p \ ~ ^Mm^, it is not appropriate to employ the frequently used heavy 
baryon formalism for the nucleons. Let us demonstrate the failure of the heavy baryon 
approach for a simple (but generic) example. Consider diagram b) in Fig.|3] which involves a 
propagating nucleon after emission of the real 7r°. We show that this nucleon propagator can 
not be expanded in powers of 1/M in the usual way. Let = (1,0) be the four- vector which 
selects the center-of-mass frame. The four-vector of the propagating nucleon is Mv + k with 
k^ = (—m n /2,p) and k 2 = —Mm n . We start on the left hand side with the correct relativistic 
result and then perform the usual 1/M expansion of the heavy baryon formalism, 



1 _ 2M _ 1 _ 1 - / -g y_ 2 

(Mv + k) 2 -M 2 vk + k 2 /2M vk^\2Mvk) m„ ^ ' ' K ' 

On sees that infinitely many terms of the 1/M-expansion contribute to the same order. The 
resulting series does not even converge and oscillates between zero and twice the correct 
answer. The source of this problem is the extreme kinematics of the reaction NN — > NNw 
with \p\ ~ y/Mm^. In that case the leading order operator = iv ■ d and the next- 
to-leading order operator = —d ■ d/2M lead to the same result, here, ±m n /2. The 
heavy baryon formalism (and therefore also HBChPT) can not cope with external momenta 
as large as \p\ ~ y/Mm n . This problem is merely related to "trivial" kinematics and can be 
immediately overcome in fully relativistic calculations. In the latter one does not presume 
that vertices and propagators can be consistently expanded in inverse powers of the large 
nucleon mass M. The occurrence of such non-expandable nucleon propagators is generic to 
iVW — > NNir and they will of course also enter in loop diagrams. In order to treat correctly 
the kinematics, one has to evaluate the loops relativistically, even though this may spoil the 
one-to-one correspondence between the loop and small momentum/quark-mass expansion. 
A consistent power counting scheme for processes iViV — > NNtc involving loops of arbitrary 
high order is not in sight at the moment. 
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We thus turn to the relativistic description of the chiral pion-nucleon system. Instead of 
giving the Lagrangian (see e.g. ref.[24]|), we mention only the relevant pion-nucleon vertices. 
The 7riViV-vertex (and also the r]NN-vertex) is of pseudovector type as required by chiral 
symmetry. The second order isoscalar chiral 7m NN contact interaction, vr°(gi) + N(p\) — > 
7r°(g 2 ) + N(p 2 ) reads, 



*2 

■J TV 



4c 1 ml+2c 3 q 1 -q 2 + — (pi+p 2 )-gi (fa +O1+P2) -?2 qi 



2M 



2M 2 



(pi+P2)-qi (pi+P2)-q2 



(9) 

This form is unique since on mass-shell it gives the most general second order (s <-> u) crossing 
symmetric polynomial contribution to the invariant isoscalar ttN amplitudes, 



A+(s,u) = ^{-4c 1 ml + c 3 (s + u-2M 2 ) + ^(s-u) 2 } , B + (s,u) 



c'2 



8M 2 



■u 



(10) 



with s and u the usual Mandelstam variables. The low-energy constants c\, c 2 , c 2 , C3 have 
already been determined (at tree level) in Jl8] from low-energy ttN data and we list their 
values for completeness: 



Cl 



-0.64, 



4 



-5.63, 



7.41 , c 3 



-3.90 



all given in GeV . As it was shown in ref. [ 25f| , the numerical values of most of these 
low-energy constants (c' 2 , c 2 , C3) can be understood largely from intermediate A excitations. 

Consider first one-pion exchange. The respective diagrams are shown in Figs. |3]a,b 
(with A4 = 7T°). We stress that these are relativistic Feynman graphs, i.e. in the inter- 
mediate states they contain the full relativistic fermion propagator which sums up several 
time-orderings. We also note that the graphs with the pion exchange after the emission of 
the 7T° from one of the proton lines does not belong to the final state-interaction according 
to our treatment (as an on-shell NN-process). We find 



^(7r,dir) 



n 3 



4M 4 (l + /i)(2 + /i) 



0.48 fm 4 



(12) 



with g n N = 13.4 the strong pion-nucleon coupling constant. We remark that it is often claimed 
that the ppn° final-state should be suppressed due to chiral symmetry. This argument is based 
on the assumption of the exchanged pion being soft, which, however, is not the case. The ir°p 
amplitude with one 7r° off its mass-shell is only of linear order in the pion mass m n (since 
k 2 = —Mm n ) and this factor of m n is cancelled by the pion propagator [m n (M + m^)] -1 . 
Consequently, as Eq.(|l^) shows, the threshold amplitude A in l7r-exchange approximation 
does not vanish in the chiral limit as it is often claimed. In fact, chiral symmetry does not 
distinguish the pp — > ppn° process from the other NN — > NNn channels. Next, there is the 
so-called n° rescattering, as shown in Fig. ||c, 



^(7r,rcs) 



9ttN H 



f 2 M(l + fx) 



£3 
2 



c + 1 + 



2ci 



0.46 fm 4 



(13) 



with /„. = 92.4 MeV the pion decay constant. Again, there is a marked difference to the heavy 
baryon case. To leading order, the relativistic couplings c' 2 , c 2 combine to give the c 2 = c' 2 + c 2 
term in the heavy baryon approach. In previous HBChPT calculations the c 2 term was found 
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with an incorrect prefactor 1/2. The relative factor of two in the relativistic calculation comes 
from the fact that products of nucleon and pion four-momenta are not dominated anymore 
by the term nucleon mass times pion energy. For the reaction NN — ► NNn the product of 
the nucleon and pion three- momenta can be equally large, since \p\ ~ y/Mm^. 

Interestingly, the combination of low-energy constants in Eq.(|l~3"D is dominated (to about 
90%) by the last term ~ — 2ci, which is related to the so-called pion-nucleon sigma-term, 
0"7r7v(O) = —4c\m\ (to leading order). Therefore the strength of the 7r°-rescattering (at 
threshold) is almost entirely due to this particular chiral symmetry breaking term. The effects 
from the A(1232)-resonance encoded in the low-energy constants c' 2 , c^', C3 turn out to be 
very small. In order to check this interpretation, we have evaluated the contributions from 
explicit A(1232)-excitations, using the Rarita-Schwinger formalism and the well-satisfied 
coupling constant relation g n ^A — ^9-kn/V^- We find: 

.(A) = 9*Nf* f M 2 (4 + /i)(M A + M(l- / i)) 

16M|M 3 (l + /i)l M| + M 2 (2/i-l) 

+(2Z - 1)[4(1 + Z)M A - M(2Z (3 + fi) + 1 + fi)] \ = 0.04fm 4 , (14) 

where we used for the off-shell parameter Z = —0.3, the value which maximizes the A(1232)- 
contribution to the P33 nN scattering volume. Such a value of Z is also consistent with neutral 
pion photoproduction off protons. The small value of A^ confirms the interpretation of the 
rescattering contribution Eq.flTB"!) given above. 

The next Goldstone boson which can contribute is the 77(547). Consider the graphs in 
Figs. |||a,b with A4 = 77. We find 

2 

A {r,) _ 9kn9 v n V - n no f m 4 (1 51 

A "4M2(m2 + MV)(2 + /^) " ' [ ' 

where we have employed the SU(3) value for the rjN coupling constant together with the 
simplified ratio of the octet axial vector coupling constants D/F = 1.5, which leads to g n ^ = 
4.6. Since this contribution is tiny, the precise value of this coupling does not matter. 



3.2 Pion loop effects 

We do not attempt a full one-loop calculation here, but rather consider certain (simple) classes 
of loop graphs which are genuine to the process under consideration. We use dimensional 
regularization and minimal subtraction to eliminate divergences and set the renormalization 
scale equal to the proton mass M. For estimating the genuine size of pion loop effect, such 
a procedure ignoring renormalization via counterterms should be sufficient and the resulting 
numbers should be considered indicative. 

Consider first a certain class of pion loop diagrams which involve the irn interaction, com- 
pare Fig. |4j. Notice that only the full class of diagrams is independent of the choice of the 
interpolating pion field. Within the calculational scheme mentioned above, this class of graphs 
gives 

A**** = #^#^(ln / ,-i + v ^lnl±^±^ 
(87rM^) 2 (l+/i) I ^ 2 V ^ 2v /7i 

1 



dx I dy— r- r , \ = -O.lOfm 4 . (16) 
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We remark that the expression in the curly brackets is not singular in the chiral limit fi — > 0. 
It has the following //-expansion: — 1/2 + ^/Jl n 2 /8 + 0(fi In fi). Again, one faces here the usual 

problem of relativistic loops. w 4.( 1o °P' 1 ) does not vanish in the chiral limit fi — > and is therefore 
not suppressed by powers of the pion mass m n compared to tree graphs. Nevertheless, its 
numerical value is small. The threshold amplitude ^4.( lo °P' 1 ) is actually proportional to the 
relativistic loop contribution |24| to the nucleon scalar form factor evaluated at an invariant 
momentum transfer squared t = —Mm n , 



j(loop,l) 



g nN {2 + fi)(fi - 1) 
12M/>2 (1 + /i) 



cr nN (— M 2 /i)ioo P 



(17) 



Whereas ref . |16| claimed that these loop graphs were sizeable and are essential for a quantita- 
tive description of the data, we find that they give only a small —4% correction compared to 
the empirical value of A. The reason for this discrepancy is to be found in the inappropriate 
application of the heavy baryon formalism to pp — ► ppn in ref. [ 16f] . Next, we consider the 
entire class of loop graphs proportional to g^N I ' fti see graphs in Fig.| and Fig.|. Figs.|| (a) 
and (b) represent -k + -k~ exchange between protons (in form of a bubble diagram) before or 
after the emission of the neutral pion. We find: 



^(loop,2) 



g 7rN (2 - /i)/i 
768vrV. 4 (2 + /i) 
-0.01 fm 4 . 



1 + 6ft) ln/t - - - 4/t + (1 + 4/i) 3/2 In 



1 + y/l + 4/i 

2^ 



Note that these loop corrections vanish in the chiral limit fi — > 0. In the class of diagrams 
proportional to g^N / ft there is also 7r°-exchange between the protons with tt + rescattering 
on one of the nucleons, see (c) and (d) in Fig.H. These graphs give 



^(loop,3) 



g-wN fJ- 



[87r/?(l+Ai)(l-2/i)] a 



(1 + /i)(l - 2/i)(4 - 6/i - fi 2 ) + (11 - 16/i + ft 2 + /i 3 ) 



2 + fi + Jfi(4 + n) 



■li In n + /t(l + fi) (2 - 2/t - /i 2 ) Jfi(A + fj) In 

= 0.22 fm 4 . (19) 

This contribution vanishes again in the chiral limit \i — > 0. Finally, there are graphs where 
a 7r + 7r - -pair is emitted from one proton and the ir° emission proceeds via charge exchange 
from the other nucleon line, see Fig.|6|. Straightforward evaluation of these diagrams gives 

g-KNli f M9/i - 2 - ii 2 ) , 6 - 12/i + 19/i 2 + 2/i 3 + // 2 

~ in jji 



^(loop,4) 



(8tt/ 2 ) 2 12(1 + /.)(! -2pt) 



2(l + /x) 2 (l-2 M )2 



/i(4 + /i) In 



2 + fj, + Jfi(A + fi) 



2(1 -2ft) 2 vr v " ' r ~ /_ " 2 
2x - 2 - 5y + /x(3 - 3x - y) + ft 2 (3 - x) 



+ 



dx dyy x 
o io 



+ 



2 - 2x - 3y + fi{x - 1) + /t 2 



-Z/ 2 + M 1 - a:)(sc - 3j/) + /i 2 (l - 2// + xy) y 2 + /i(l - a?) (a; + y) + fi 2 (l - y) 

= -0.25 fm 4 , (20) 

a contribution which again vanishes in the chiral limit. Summing all these various loop 
contributions in Eqs . flIB| 18| , [19| , p0D , we find 

^ (loop) = -0.14 fm 4 , (21) 

which is a rather small number as compared to the tree level pion exchange or the vector meson 
contributions (section 3.3). We observe sizeable cancelations between the various classes of 
loop graphs and conclude therefore that these do not play a significant role for explaining the 
threshold amplitude A. The imaginary part of A is discussed in appendix B. 
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3.3 Heavy meson exchanges 

Due to the large momenta involved, one expects additional contributions due to the exchange 
of heavier mesons, which are of much shorter range than the pion exchange considered so far. 
From the symmetry point of view such terms are much less constrained and thus exhibit a 
certain unavoidable model-dependence. Nevertheless, these terms can play a significant role 
as first stressed by Lee and Riska J7]] and confirmed by Horowitz et al. ||. We do not discuss 
here additional effects from vector meson nucleon form factors, which are a model-dependent 
(and unobservable) concept to account for the finite size of the hadrons involved. Note also 
that in a strict field-theoretical sense such form factors can not be uniquely defined. At the 
invariant momentum transfer t = —0.127 GeV 2 the form factor effect is not expected to be 
large. For example, one expects a 10% effect for typical monopole form factors with cut-offs 
~ 1.5 GeV. 

Consider first neutral vector mesons. We start with the a; (782). There is sizeable un- 
certainty about its coupling constant to the nucleon, extreme values are e.g. found in the 



dispersion-theoretical analysis of the nucleons electromagnetic form factors, g^N — 21 127 . 
However, it is not clear how one has to transcribe such a value to the one-boson exchange pic- 
ture of the NN force. In conventional boson-exchange models, the inclusion of the correlated 
up continuum allows one to work with a coupling constant that is compatible with the SU(3) 
value, g^N — 9 or the value g^N = 10.1 ± 0.9 found from forward NN-dispersion relations 



f28| . For a detailed discussion, see e.g. ref. [29|. There is agreement that the tensor-to-vector 



coupling ratio of the cu-meson is very small. If we set k w = and use the coupling constant 
9wN = 10, we get 

A {») = 9nN9l N (2-fj,) = 4 

2M 2 (M 2 + M 2 p)(2 + p) ' V ; 

which is quite sizeable. In a similar fashion, we evaluate the p(770) contribution. Here, 
there is less debate about the coupling constant g p ^ and also, it is well established that the 
tensor-to-vector coupling ratio k p is large. For simplicity, we use g P N = 3 (obtained from 
p-universality, g P N = g p /2 with g p = 6) and k p = 6. That leads to 

AM = 2M 2 (M/rA)(2 + ,) [ 2 + ^ - «> - 1) + ^ + I*)] = ° ^ • < 23 > 

We note that this form is very different from what has been used in the literature so far, where 
one finds the u, p-exchange to be proportional to /i(l + k^^). We do not employ any inappro- 
priate non-relativistic approximation here and thus obtain the results shown in Eqs.( p2] , [2"3| ). 
Note that due to the large tensor-to-vector coupling ratio of the p-meson k p = 6 the terms 
proportional to p and p 2 in the square bracket of Eq.(|23|) are most important. We have also in- 
vestigated the role of 0(1020) exchange, which can be related to the strangeness content of the 
nucleon wave function. The values of the <pN coupling constant span a large range, as docu- 
mented in table 2 of ref. |30| . To get a more precise number, we proceed as follows. We assume 
that the tensor-to-vector coupling ratios are given by the dispersive analysis of the nucleon 
electromagnetic form factors, k p = 6.1, k u = —0.16 and = —0.22 f2"7|| . Using furthermore 



the SU(3) relation g^N = 3gpN — V2g<f>N 123 together with g pN = 2.63 |31[] and g^N — 10-1 p8 
we have g+ N = -1.56. This leads to A {Y) = +A {uj) +A W = (0.38 + 1.40 + 0.02) fm 4 = 1.80 
fm 4 , which is practically identical to the result obtained above using the simplified coupling 
constants. Thus, the sum of the vector meson contributions is fairly stable against parameter 
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variations and also the 0(1020) does not play any role. Finally, we remark that substituting 
(g P N, M p , k p ) by (e, 0, n p ) one can convince oneself that the one-photon exchange gives a tiny 
correction of A^ = 0.0085 fm 4 , as it is of course expected. 

Another mechanism, first proposed in this context in fli~3"| |, is the emission of the neutral 



pion from the anomalous uipn- vertex, with the p(770) coupling to one and the a; (782) to 
the other proton. This type of graph is similar to well-known meson-exchange currents in 
electromagnetic processes. The pertinent interaction vertex, with its strength given by the 
coupling constant G^ p7r , can be determined from the anomalous Wess-Zumino-Witten term 
for vector mesons (we only show the part of relevance here), 

= <r# {dp Uv ) p a ■ d P 7? , (24) 

with e^ va P the totally antisymmetric tensor in four dimensions (e 0123 = — 1). Using again the 
universal p-coupling, g p = 6, the gauged Wess-Zumino-Witten term for vector mesons leads 
to the coupling constant 

Gujpn = = 1-37 • (25) 

87T Z 



Similar values for G upiT have been found in refs. p2|, [33|j from systematic studies of u;(782) 



and 0(lO2O)-decays. The contribution of the anomalous cu/wr-vertex to A is given by 
A ( Uf m) _ 9u J n9 p n{1 + ^G^Mp 2 4 

A ~ 2UMI + M»(M 2 + M-p) ~ °- U9im ' {2b) 

which is at first glance quite small. However, due to the two derivatives appearing in Eq. (E^) , 
one expects this contribution to be of much larger importance in the P-wave amplitudes. 
An analogous short range mechanism is the ir° emission from the ao^-vertex, proposed in 
ref.[|n|. We have evaluated the respective contribution to A and found that it is negligibly 



small, A^ aor! ^ = 0.003 fm 4 . Finally, we do not consider a scalar meson exchange here. First, 
there is no scalar meson resonance which couples strongly to the nucleon and secondly, the 
fictitious "<j(550)" of one-boson exchange models just simulates the long and intermediate 
range part of uncorrelated 27r-exchange in the NN interaction ||34|| . The latter comes along 
with pion loops and is to some extent contained in the loop graphs considered above. 



3.4 Total threshold amplitude 

We are now in the position to evaluate the full amplitude A from the various contributions. 
Combining Eqs. (0,0,0,0,^,^,^) , we get 

^(thy) _ ^(7r,dir) _|_ ^(yr.res) _|_ ^(77) _|_ ^(loop) _|_ jj^ui) _|_ j^p) _|_ J^pir) 

= (0.48 + 0.46 + 0.02 - 0.14 + 1.35 + 0.48 + 0.09) fm 4 

= 2.74 fm 4 , (27) 

which compares well with the empirical value given in Eq.(Rl). The resulting total cross section 
Ctot(Tiab) is shown in Fig. || by the dashed-dotted line. Of course, taken the uncertainty in 
certain coupling constants and our simplified treatment of the final-state interaction, the 1% 
agreement between our theoretical prediction for A and its empirical value, Eq.(|7]), should 
not be taken too serious. We only want to make the point that these well-known boson- 
exchange diagrams, when evaluated fully relativistically, can explain the near threshold data 
for pp — > ppir . Note also that the contribution from explicit A(1232)-excitation, cf. Eq.(|i~4T), 
is contained in the it rescattering term via the low-energy constants c' 2 , c^, C3. 
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4 Charged pion production in ^—collisions 



In this section we will discuss charged pion production using the same approach. The T- 
matrix for charged pion production in proton-proton collisions, pi(p )+P2{~P ) ~~ > p+n+ir + , at 
threshold in the center-of-mass frame reads (see also appendix A for the general NN — > NNn 
threshold T-matrix), 

T?h (pp -> pnn + ) = ^(iff 1 -i? 2 + ff 1 xff 2 ). p - V2Bi(a x + a 2 ) ■ p . (28) 
V 2 

As indicated by the ordering of the three particles pnir + in the final state, the spin-operator d\ 
is understood to be sandwiched between the spin-states of the ingoing proton pi(p) and those 
of the outgoing proton, while <?2 acts between the proton P2{—p) and the outgoing neutron. 
Since the pn-system at rest can be in a spin-singlet or in a spin-triplet state two different 
transitions are possible at threshold. The threshold amplitude for the singlet transition 3 P — > 
1 Sqs is by isospin symmetry proportional (with a factor 1/V2) to the threshold amplitude 
A for the reaction pp — > ppir introduced in Eq.(|I]). The new threshold amplitude for the 
triplet transition 3 Pi — > 3 Sis is called B and the factor — v2 was taken out for convenience 
in Eq . (|28|) . In analogy to Eq.(3) one deduces from unitarity 

B = \B\ e l5(3pi) , (29) 

with the 3 Pi pp phase shift to be taken at T^ h = 292.3 MeV, where 5( 3 Pi) = -28.1° (FA95 
solution of VPI). Because of this larger phase fixed by unitarity the imaginary part ImB 
will contribute non-negligibly to the total cross sections near threshold. Notice, that the 
abovementioned unitarity relation Eq.(|29"D neglects the (small) inelasticity due to the pp — > 
dn + channel which opens 4.8 MeV lower at Ti a b = 287.5 MeV. 



4.1 Extraction of the threshold amplitudes 

Employing the same method as in section 2.2 to correct for the strong S-wave pn final-state 
interaction, the unpolarized total cross section for pp — > pmr + reads, 



/Mx3 2^ 

0tot(-7lab) = ' ' 



V4tt^ (2M + T lab )3/2 
x / liTV A (^ 2 ' M2 ' M n) A (^ 2 ' rnl + AM 2 + 2MT lab ) 

JM+M n W v 

{\A\ 2 F s (W) + 2\B\ 2 F t (W)} . (30) 



x 



The correction factors from the pn singlet and triplet S-wave final state interaction are given 
in the effective range approximation by 

F s ,t(W) = {l + a s , t (a s ,t + r s , t )Pl + \al t rl t Pt] ' , (31) 



with W the final-state proton-neutron invariant mass and W max — V4M 2 + 2MTi ab — . 
M still denotes the proton mass and M n = 939.57 MeV stands for the neutron mass. The 
quantity P 2 = \(W 2 , M 2 , M 2 ) / 4W 2 is the squared pn center-of-mass momentum. The singlet 
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and triplet scattering lengths and effective range parameters for elastic rip-scattering are taken 
from g2j, their empirical values being a s = (23.748 ± 0.010) fm, a t = (-5.424 ± 0.004) fm, 
r s = (2.75 ± 0.05) fm and r t = (1.759 ± 0.005) fm. We neglect here the coupling between the 
3 Si and the 3 Di pn-states, which should be very small at the energies under consideration. 
Such effects go beyond the accuracy of the effective range approximation and our treatment 
of the pn final-state interaction. 

The data base of total cross sections for the process pp —>■ pmr + in the 30 MeV region 
above threshold consists at present of five data points measured at IUCF @]. We leave out 
the data point at the highest energy T lab = 319.2 MeV where the n + angular distributions are 
no more isotropic and thus P- waves start to become important. We also found it important 
to ignore the data point at the lowest energy Ti a b = 294.3 MeV. Using Eq. (|3~0|) for the total 
cross section and the value of |^4| = 2.72 fm 4 as determined from the pp — > pptr data, one 
finds in a least square fit of the remaining three data points for the modulus of the triplet 
threshold amplitude, 

|£|=3.16fm 4 , (32) 

with a very small total \ 2 = 0.044. For these values of |^4| and \B\ the data point at 
^lab — 319.2 MeV is underestimated by 17% in the S-wave approximation (see table 2). An 
unconstrained fit of the same three data points gives |^4| = 3.00 fm 4 and \B\ = 3.15 fm 4 with 
a marginally smaller total \ 2 = 0.042. It is quite remarkable that |^4| is found to be in 10% 
agreement with the value obtained from fitting the many precise near threshold pp — > ppn 
data. Due to the very strong pn final-state interaction in the 1 Sq exit channel (a s = 23.75 
fm) the singlet transition contributes a factor 30 to 40 less to the total cross section than the 
triplet transition. From our fit we get at the lowest energy Ti a b = 294.3 MeV a total cross 
section of 0.57/zb, compared to the experimental value of (0.71 ± 0.04) fib given in ref. Q. 
Only if one widens the error band of this data point by the ±15% absolute normalization 
uncertainty it becomes (at the lower end) marginally consistent with the remaining data. If, 
however, the data point at Ti a b = 294.3 MeV were included in an unconstrained fit, very 
different values of the threshold amplitudes, |^4| = 6.76 fm 4 and \B\ = 2.91 fm 4 , would be 
found with a total \ 2 = 1-0- I n particular |^4| would be a factor 2.5 larger than the one 
obtained from the fit to the pp — > ppn data. Of course, such a large deviation from isospin 
symmetry is unacceptable. 



T !ab [MeV] 


294.3 295.1* 298.0* 299.3 306.3 314.1 319.2 




0.71 1.1* 3.84* 4.81 13.91 25.5 41.1 
±0.04 ±0.3* ±0.37* ±0.24 ±0.65 ±1.6 ±1.7 


w 


0.57 1.07 3.46 4.81 13.82 25.76 34.14 



Table 2: Total cross sections for pp — > pnn + as a function of 71 a b. The data are taken from 
ref.[|J and the preliminary data points from COSY-TOF [35] are marked by an asterik. The 
(constrained) fit is described in the text. 



Anticipating the positive sign of the real part from the calculation in the following section 
and using the information from the 3 Pi phase shift, one gets the following experimental value 
of the triplet threshold amplitude B, 

B (ex P ) = ( 2 .8-zl.5)fm 4 . (33) 
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This number should be considered indicative since the systematic error of the extraction 
method is not under control when only three data points are fitted. Notice also that the 
imaginary part is fairly sizeable, quite different from A, the threshold amplitude for pp — > 
ppir . 



4.2 Diagrammatic approach 

Next, we turn to the evaluation of the relativistic Feynman diagrams contributing to pp — > 
pnii + at threshold. In addition to the ones considered for pp — > ppir , there is now the 
possibility for isovector pion-rescattering. The chiral ttN Lagrangian contains such vertices 
at leading order (the so-called Weinberg-Tomozawa vertex) and at next-to-leading order (a 
vertex proportional to the low-energy constant C4 = 2.25 GeV^ 1 ). In fact there are always two 
isovector rescattering diagrams, one with a 7r° and another one with a tt + being exchanged 
between the nucleons. The isospin factors of both diagrams are equal with opposite sign. 
However, before adding them one has to account for the fact that the role of p and n is 
interchanged in both graphs. This is done by multiplying those graphs where the final state 
neutron n comes from the initial state proton pi(p) with the negative spin-exchange operator 
— (1 + d\ ■ <?2)/2. Altogether, one finds from isovector pion-rescattering at leading and next- 
to-leading order, 

B( .,^ w^i) = _ . 82tmV (34) 

Of course, we neglect here the small isospin breaking due the different charged and neutral 
pion masses and the different proton and neutron masses. ^From the other pseudoscalar meson 
(tt and rj) exchange diagrams, one finds the following contributions to the triplet amplitude, 

w(7r,dir) _ ^iv( 3 + 2 / i ) -I58f m 4 (W) 

£(-,res) = ^(^s) = .46fm 4 , (36) 
B (v) = _ A (v) = _ .02fm 4 , (37) 

(see also Eqs. ([13|JT5|) ). Note that the (direct) l7r-exchange contribution B^' dl ^ is rather large 
due to an enhancement factor 3+2/x in comparison to given in Eq. (|T2"| ) . The dominance 

of the chiral one-pion exchange in the triplet transition 3 Pi — > 3 S\s found here supports the 



argument of [15| concerning the dominant role of the long range (chiral) pion exchange in the 
reaction pp dir + which proceeds via the same transition near threshold. 

^From the vector meson (p and u) exchange diagrams one finds the following contributions 
to B, 

B (ul) = TT^TT^^fr? r = 1.46fm 4 , (38) 

M 2 (M2 + M 2 /i)(2 + /i) ' v ; 

= 9jNgJ N (^ P - 4) r , ^ = _ 
4M 2 (M p 2 + M 2 /i)(2 + /i) L P V 4/J 



B {wpn) = . (40) 

One observes that the sizeable tu-exchange contribution B^' is approximately equal to A^ 
given in Eq.(E^). Note also that the expression for tu-exchange cannot be recovered by simply 
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substituting the vector meson mass and coupling constants in the expression for p-exchange. 
The reason for that are certain diagrams with (charged) p + -exchange which have no analogy 
in the case of the (neutral) w-meson. In the spirit of vector meson dominance one could also 
think of an isovector npNN contact vertex of the form 

ArpiV = ^N^r • (vr x p, )N . (41) 

The form of this vertex and the coupling constant in front are copied from the Kroll-Ruderman 
term for charged pion photoproduction replacing the charge e by the universal p-coupling 
g p = 6 and the photon field by the isotriplet p-meson field p^. The respective p° and p + 
exchange diagrams give rise to the following contribution to the triplet amplitude, 

B (p,KR) = 9-N9 P 9 P N _ M \ = 45f 4 (42) 

M 2 (M 2 + M 2 p)V 4 p ) y 1 

Admittedly, this contribution is somewhat speculative since the Kroll-Ruderman vertex for 
p-mesons presumes a particular realization of vector meson dominance. We do not investi- 
gate in further detail pion loop diagrams contributing to B, since they will turn out to be 
small in analogy to the amplitude A discussed in section 3.2 (for the classes of diagrams 
calculated there). The imaginary part of the triplet amplitude B can only be generated in 
one-pion loop approximation by the two-pion exchange box diagrams shown in Fig. 9 and 
evaluated in appendix B. It can, however, not be expected that the one-pion loop approx- 
imation will be sufficiently accurate, since the process NN — > NNtt at threshold is also 
sensitive to short distance dynamics. This needs further study but goes beyond the scope 
of this paper. A detailed discussion of the imaginary parts of A and B in one pion loop 
approximation is given in appendix B. Summing up the various tree level contributions given 
in Eqs.(H^5l^0|^^|gPD we get, 



B (thy) = 2.74 f m 4 , (43) 

which is very close the real part of the experimental value in Eq.(^), Rei3 ( - exp ^ = 2.8 fm 4 . We 
thus conclude that also the real part of the triplet amplitude Re B can be well understood in 
terms of these well-known tree-level meson exchange diagram when evaluated fully relativis- 
tically. Of course, the relatively large imaginary part Imi3^ cxp - ) = —1.5 fm 4 (more than half as 
large as the empirical real part) remains unexplained in tree approximation. As mentioned 
earlier this imaginary part originates (because of unitarity) from the fact that the 3 Pi pp 
phase shift is rather large at the pion production threshold. The large imaginary part ImB 
thus reflects the strong initial state interaction in the 3 Pi entrance channel. For the singlet 
transition amplitude A the situation is different. Accidentally, the 3 Po pp phase shift is very 
small at the pion production threshold and thus there is only weak initial state interaction in 
the 3 Pq entrance channel. Another mechanism which could contribute (significantly) to Im£> 
is the two-step process pp — > dir + — > pmr + . The threshold for the deuteron channel pp — > dir + 
opens 4.8 MeV lower at Ti a b = 287.5 MeV and the corresponding total cross sections are more 
than an order of magnitude larger [[| || than the ones for pp — > pwK + . 

We have done several checks on our treatment of the pn final-state interaction. First, we 
compared the effective range approximation with the empirical values of the 1 5'o and 3 5*i pn 
phase shifts in the energy range relevant here and found that deviations are smaller than 2%. 
Secondly, we have studied the tt + production cross section as a function of the pn invariant 
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mass W (given by Eq. (|30"D without rfVT-mtegration) and found good agreement with the data 
of 0j , cf . Fig.[?|. ^From all this we conclude that our treatment of the pn final-state interaction 
is fairly realistic. 

Finally, we believe that all the features of the processes pp — > pNn that we have learned 
here from the relativistic diagrammatic approach presented here will be useful for further 
studies within more complete dynamical models, which e.g. treat initial- and final-state in- 
teractions simultaneously. We also remark that a similar covariant one-boson exchange model 
has been developed in ref. J36|, which describes the data at much higher energies, 0.3 GeV 



< TJab < 2.0 GeV. The same model has been applied to threshold data in ref.|p7||. These au- 
thors introduce in addition (unobservable) meson nucleon form factors and energy-dependent 
coupling constants. This makes a direct comparison between their work and ours very diffi- 
cult. However, no isoscalar pion rescattering and no pion loops are considered. Similar to our 
finding it is concluded that u;(782) exchange is important close to pion production threshold 
and that the A(1232) plays no role. 



5 Eta— meson production in ^—collisions 

In this section we will discuss //-production using a similar approach. The T-matrix for 
77-production in proton-proton collisions, Pi{p) + P2{—f>) ~~ > P + V + Vi a ^ threshold in the 
center-of-mass frame reads (see also appendix A), 

T th (PP PPV) = C (i o x - i a 2 + o x x a 2 ) ■ p. (44) 

with C the (complex) threshold amplitude for //-production. The //-production threshold is 
reached at a proton laboratory kinetic energy T x % = m v (2 + m v /2M) = 1254.6 MeV, where 
m„ = 547.45 MeV denotes the eta-meson mass. 



5.1 Extraction of the threshold amplitude 

In the case of //-production near threshold it is also important to take into account the r]p 
final-state interaction, since the //iV-system interacts rather strongly near threshold. In fact 
a recent coupled-channel analysis ]38| of the (7rN, r)N)-system finds for the real part of the 



r/N scattering length Rea^ = (0.717 ± 0.030) fm. For comparison, this value is a factor 5.7 
larger than the n~p scattering length, a n ~ p = 0.125 fm [jS^I , measured in pionic hydrogen. 

We assume that the correction due to the S-wave rjp final-state interaction near thresh- 
old can be treated in effective range approximation analogous to the S-wave pp final-state 
interaction. We furthermore make the assumption that the final state interactions in the pp 
subsystem and in the two r\p subsystems do not influence each other and that they f actor ize. 
These are of course very strong assumptions, but as we will see soon, such a simple ansatz for 
the final-state interaction in the ppr] three-body system allows to describe rather accurately 
the energy dependence of the total cross section a to t{pp — ► ppr]) near threshold. Using the 
factorization ansatz mentioned before, the unpolarized total cross section for pp — > ppr] reads, 



x f Sv dsrj F v (s v ) F V (6M 2 + 2MT lab + m\ - W 2 - s v ) . (45) 
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with Wmax — V4M 2 + 2MTi a b — the endpoint of the di-proton invariant mass spectrum 
and Fp{W) given by Eq.(6). The variable s v is the invariant mass squared of the first 7?p-pair 
and the argument of the last function F v (s v ) in Eq. fl45|) , s v = QM 2 + 2MT lab + m 2 — W 2 — s v , 
is the invariant mass squared of the second ??p-pair. The expressions 



3M 2 + MT lah + -(m 2 l -W 2 )± 



1 



2W 



{W 2 - 4M 2 )X{W 2 , m 2 4M 2 + 2MT lab ) (46) 



give the boundaries of the ppr\ three-body phase space in the (s v , iy 2 )-plane. Obviously, 
the formula for the total cross section, Eq. 



s v , since 



is invariant under the permutation of the 

due to the 



two rjp-paiis, s v 

S-wave 7]p final-state interaction reads in effective range approximation, 



Furthermore, the correction factor F^s^ 



Frj (Sjj) 



/r)iv( S ??) 


2 







1 



^\( Sr] ,m 2 ,M 2 ) + ^p^X(s v , m% M 2 ) 



8s r 



(47) 



Here, f^( s v) * s ^ ne S-wave r]N elastic scattering amplitude. The (complex) r/N scattering 
length a vN = ((0.717 ± 0.030) + i (0.263 ± 0.025)) fm is taken from ref.0 and the (complex) 
t]N effective range parameter r v ^ = ((— 1.50±0.13) — i (0.24±0.04)) fm stems from ref . ||40|| . It 
is important to note that both ref. |38j and ref. |4(J using quite different methods agree within 
error bars on the value of the r] N scattering length a V N . 

Using Eq. ( ^5| , ^6| , ^7D for the total cross section and the central values of a n ^ and r vN one 
finds in a least square fit of the six data points from CELSIUS [[J for the modulus of the 
threshold amplitude 

|C| = 1.32fm 4 , (48) 

with a total x 2 = 3.8. The resulting energy dependent cross section from threshold up to 
^lab = 1375 MeV is shown in Fig. |] together with the data from CELSIUS p^.p 2 It is rather 
astonishing that one can describe the total cross section data up to 100 MeV above threshold 
with a constant threshold amplitude C and a simple factorization ansatz for the three-body 
final-state interaction. 



5.2 Diagrammatic approach 

Next, we turn to the evaluation of the relativistic Feynman diagrams contributing to pp — > 
at threshold. The resulting expressions can essentially be copied from the case pp — > i 
making only the substitution (g^Ni™-*) — > (grjN, m ri)- One finds for n°, t], u>, p°-exchange 



£(»7,dir) 



9vN9iTN m V 



4M 2 (m 2 + Mm v )(2M + m v ) 

9^N 

4M 2 (M + m ri )(2M + m v ) 
9vn91n( 2M ~ m v) 



2M 2 (M 2 + Mmr,)(2M + m, 



0.17fm 4 
0.02 fm 4 , 
- = 0.23 fm 4 



9yNg pN 



2M(M 2 + Mm v )(2M + m v ) 



x 



2 + 



1) + ^( 1 + 8 K ' 



(49) 
(50) 
(51) 



0.50 fm 4 



(52) 



#2 



Earlier data from SATURNE ||, ||, || do not have the same accuracy and are not considered further. 
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Note that the p° exchange has become dominant because of the large tensor-to-vector cou- 
pling ratio k p = 6 and the larger ratio m v /M = 0.58.0 Besides these diagrams with rj- 
emission before and after meson exchange between the protons, one has to account for the 
strong rjp rescattering. Microscopically, the strong rjN S-wave interaction originates (among 
other things) from the nucleon resonance Sii(1535) which is supposed to have a very large 
coupling to the 7/iV-cliannel. Instead of introducing this resonance together with several 
parameters (mass, width, coupling constant), we will merely introduce here a local NNr/rj 
contact vertex of the form 

C vN = KN(x)N(x) 7] 2 (x) . (53) 

The interaction strength K is then determined by the real part of the rjN scattering length. 
This means that the pseudovector Born graphs plus the contact vertex sum up to give the 
empirical value of Rea^ = 0.717 fm. This leads to the equation 

4tt(i + ^f) Rea vN = 2K - „ r ,^f\^ . (54) 



M I r,JV M(4M 2 - mj) 

which results in a value of K — 7.41 fm. The ^-rescattering graph (analogous to Fig. ^c) 
leads to the following contribution to the threshold amplitude, 

C („,res) = 9jjnK_ = Q 4Q 

Evidently, all contributions to C scale with the (empirically not well determined) r/iV-coupling 
constant g^. The numbers given in Eqs. (45,46,47,48,51) which add up to the empirical value 
of \C\ = 1.32 fm 4 follow with g V N = 5.3. Such a value of g V N is consistent with all existing 
empirical information on it. The SU(3) flavor symmetry connects the pion-nucleon and eta- 
nucleon coupling constants via the D/F ratio (of the baryon octet axial vector couplings), 

3 -D/F , , 

^ = 9 * N V3(1 + D/F) ■ (56) 

Using g nN = 13.4, the for our purpose optimal value g v ^ = 5.3 requires a ratio D/F = 1.37. In 



fact a systematic analysis of semileptonic hyperon decays in ref. |45j gives D/F = 1.58 ±0.07, 
not far from this number. Of course, SU(3) is broken to some extent by the strange quark 
mass. For a recent update, see e.g. jRJ. There is one further contribution we have not 
discussed so far, related to the ao?7vr-coupling. This coupling is rather uncertain. If one 
assumes that the diagrams with ^-emission from the a ?]7r-vertex contribute with a positive 
sign, 

2 r> 9\ 



c (aoipr) = 9a NgnN(c d m, n - 2c m „^ ; = 05 fm 4 (57) 

x/6M/2(M a 2 o +Mm v )(ml + Mm v ) ' 

one can even lower the r/A^-coupling constant to g V N = 5.1. We used here g ao N = 4.5 
and Cd = 32MeV, c m = 42MeV j48|. We also note the one-boson exchange model of 
for elastic NN-scattering uses an r^iV-coupling constant of g9^ B = 6.8 not far from our value 
Styv = 5.3. The results of our approach applied to the reaction pn — > pnrj are briefly discussed 
in appendix A. 

The main point we want to make here is that even the pp — > ppr] threshold amplitude 
can be understood in terms of these well-known meson exchange diagrams when evaluated 
relativistically. With rather mild assumptions on the coupling constant g V N and the form of 
the 7/iV-rescattering one can easily reproduce the empirical value \C\ = 1.32 fm 4 . For another 
boson-exchange approach to 77-production emphazising the role of the Sii(1535) nucleon 
resonance, see e.g. jID], |50| (and references therein). 

# 3 Interestingly, recent measurements of the angular distributions in pp — > pprj suggest the dominance of 
vector meson exchange H . 



ml 



47 



47 
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5.3 Comments on rj'— production near threshold 



Finally, we like to comment on the recent ^'-production data near threshold from COSY flSl 



and SATURNE |43| . Taken face value, the energy dependence of the four COSY cross section 
data points is best described by the pure three-body phase space behavior (as shown by the 
dotted curve in Fig. 2 of ref.|51]]). Of course, it is hard to imagine that the pp final-state 
interaction does not play a role for pp — > pprj' that close to threshold where these data were 
measured. We have analyzed the combined COSY and SATURNE data (six data points) 
within our approach including the pp final-state interaction in effective range approximation. 
Only if one ignores the COSY data at the two lowest energies (at cm excess energies of 1.5 
and 1.7 MeV, respectively), one can fit the remaining four data points with a modulus of the 
threshold amplitude of \C \ = 0.21 fm 4 with a small total y 2 = 2.4. This fit leads to values of 
the total cross section at the two lowest energy points of 5.2 nb and 6.3 nb (see also table §). 
These numbers are more than twice as large as the corresponding central values given in 



ref.[[51[]. Note, however, that in Fig. 2 of ref.[[51[] sizeable uncertainties in the excess energy are 
given, e.g. the lowest point is at Q = (1.5 ± 0.5) MeV. Our discussion always refers to the 
central Q values. 



Q [MeV] 


1.5 


1.7 


2.9 


3.7 


4.1 


8.3 


&tot [nb] 


2.5 ±0.5 


2.9 ± 1.1 


12.7 ±3.2 


19.2 ±2.7 


25.2 ±3.6 


43.6 ±6.5 


°& H 


5.23 


6.29 


13.34 


18.45 


21.09 


50.13 



Table 3: Total cross sections for pp — > ppv( as a function of the cm excess energy Q = 
\/AM 2 ± 2MTi a b — 2M — M v > (only the mean value is given). The data are from refs. f43|, |5T| . 
The fit is described in the text. 



Within the relativistic one-boson (jr°, i], u, p°) exchange model the fit value \C \ = 0.21 fm 4 
implies the relation ^^(1 — 1.28 e) = 1.12. Here g v >N denotes the 7/iV-coupling constant and 
e is the fraction of pseudoscalar ?/iV7V-coupling (since the r/(958) is no Goldstone boson there 
is no reason to favor the pseudovector coupling). Interestingly, only the tensor interaction of 
the p-exchange (~ n p ) is sensitive to the parameter e. 

According to ref. [513] one can relate the r/iV-coupling constant to the quark helicity con- 



tribution to the proton spin, AS = y3/2g 1] i N f 7r /M + 0.15. Using the recent determination 

AS = 0.45 ± 0.09 of ref. [53j from deep inelastic polarized lepton scattering, one can extract 
within the relativistic one-boson exchange model an 7/iV-coupling constant of g^N = 2.5±0.7 
and a pseudoscalar coupling fraction of e = 0.4 ± 0.1. It remains to be seen whether other 
r/-production processes (e.g. photoproduction jp — > rj'p) are consistent with these values. 
The exchange of tc°, uj and p° contribute to \C \ = 0.21 fm 4 approximately 30%, 20% and 50%, 
respectively. For another l7r-exchange model to ^'-production, see e.g. 
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A General threshold T— matrices 



In this appendix, we write down the general form of the threshold T-matrices for NN — ► NNir 
and NN — > NNr] using the isospin formalism for the two-nucleon system. In the case of 
(isovector) pion production one has, 

T™(NN^NNtt) = ■^( l B 1 -ia 2 + a 1 x a 2 ) ■ p (n + r 2 ) ■ x* 

+ | {ai + a 2 ) ■ p (i n - i f 2 + n x f 2 ) ■ x* ■ (A.l) 

with x the three-component isospin wave function of the final state pion, e.g. x = (0> 0, 1) f° r 
7T° production and x — (l>2,0)/\/2 for 7r + -production. The complex amplitude A belongs 
to the (singlet) transition 3 P — > 1 S s with conserved total isospin 7 tot = 1. Similarly, the 
complex amplitude B belongs to the (triplet) transition 3 P X — > 3 5is, also with total isospin 
J t ot — 1- In fact the selection rules which follow from the conservation of parity, angular 
momentum and isospin allow only for these two transitions. Note also that there is an 
invariance under the substitution (A, ai,a 2 ,p) <-> (B,ti,t 2 ,x*). 
In the case of (isoscalar) eta production one has, 

T C ™(NN^NN V ) = -{ia 1 -ia 2 + a 1 x a 2 ) ■ p (3 + n • r 2 ) 

D 

+ — (i o\ % <7 2 G\ x <j 2 ) • p (1 - fi • r 2 ) . (A.2) 

Again, the complex amplitude C belongs to the (singlet) transition 3 P — ► 1 SoS with total 
isospin /tot = 1 and the complex amplitude T> belongs to the (triplet) transition 1 P\ — > 3 S\s 
with total isospin / tot = 0. The determination of the latter amplitude T> requires measure- 
ments of the total cross sections for the process pn — > pni] having neutrons either in the target 
or the beam. Finally, we note that the expressions for T^(NN -> NNn) and T^(NN -> 
NNrj) incorporate the Pauli exclusion principle, since combined left-multiplication with the 
spin-exchange operator (1 + o\ • a 2 )/2 and the isospin-exchange operator (1 + f\ • f 2 )/2 re- 
produces them up to an important minus sign. 

For the sake of completeness, we give also the contributions of the one-boson exchange dia- 
grams discussed in section 5.2 to the (triplet) 77-production amplitude V. From pseudoscalar 
meson (ir, rf) exchange one finds, 

pto = -3 C ( " 0) = -0.52 fm 4 , = C (ri ' dir) + C (ri ' rcs) = 0.42 fm 4 , (A.3) 

(see Eqs. (49,50,55)). In addition vector meson {p,u) exchange gives rise to the terms, 

(w) = 2M + m v M = 

2M-m„ ' K J 



V {p) = 



3g V N9pN 



2M(M 2 + Mm v )(2M + m v ) 



p 

x 



-2 + ^ + k,-1) + ] ^k,(i + -k, 



1.50fm 4 , (A.5) 
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(see Eq.(51)). These contributions sum up to = 1.81 fm 4 . The expression for the 

total cross section <r tot (T lab ) of the reaction pn — > pnn is obtained from Eq.(45) replacing the 
factor 2\C\ 2 F p (W) by \C\ 2 F S (W) + \V\ 2 F t (W) (see Eq.(31)) and we have neglected the small 
neutron-proton mass difference. 

Recently, total cross sections for pn — > pnn have been extracted from measurements of the 
process pd — > ppnn at CELSIUS A fit of the two data points closest to threshold (at 

Ti a b = 1296 and 1322 MeV) gives for the empirical triplet //-production amplitude |P|( cxp ) = 
2.3 fm 4 . Compared to this value the abovementioned theoretical prediction, J)^ thy ' = 1.8 fm 4 , 
is only about 20% too small. 



B Imaginary parts from one— pion loop graphs 



In this appendix we will give explicit expressions for the imaginary parts of the pion pro- 
duction threshold amplitudes Im^4 and ImS as they arise in one-pion loop approximation. 
To that order any non-vanishing imaginary part can only come from those one-loop dia- 
grams which involve proton rescattering in the initial state. These are just the two-pion 
exchange box diagrams shown in Fig. 9. We apply the Cutkosky cutting rules to evaluate 
their imaginary part. It is then given by the product of the upper tree-level subgraph (i.e. the 
threshold pion production amplitude with l7r-exchange including the leading order isovector 
pion rescattering) and of one-half the invariant pp two-body phase space times the lower tree- 
level subgraph. The latter two factors combine to the 3 Po or 3 P X pp phase shift calculated 
perturbatively in l7r-exchange approximation at = m n (2 + m 7T /2M). Altogether one gets 
thus from perturbative unitarity, 



lmA = A^ dir) -5 ln ( 3 P ) , 



(B.l) 
(B.2) 



with A {n ' dir \ and £^< dir ) given in Eqs.(|TJ| 



5p. Using the projection formulas of 
ref.[Q] to calculate the l7r-exchange 3 Po,i phase shifts one finds the following analytical 
results, 



Im.A 



9*n J + 
647rM 4 (l + y u)(2 + /i) i 



, ln ( 2 + -) 



0.54 fm 4 



(B.3) 



ImB 



<WM4 + ^) 



64ttM 4 (1 - 
M -4 



x 



/i)(2 + /i) a 



2(4 + f i) (4 + /x) 2 



3 + 2/i 



2^ 2 



In 



2 + i 



(2 + // 



-0.24 fm 4 



(B.4) 



with qa = gnNfn/M. One observes that in one-pion loop approximation Im^4 is too large 
by a factor of 2 with the wrong (positive) sign, whereas ImB is to small by a factor of 6. 
Clearly, there is important short range NN-dynamics missing in one-pion loop approximation 
as can be seen by comparing the empirical 3 Po,i phase shifts at = 77^(2 + m 7r /2M), 
<5( 3 P ) = -6.3° and 5( 3 Pi) = -28.1°, with the lvr-exchange approximation, <5^( 3 P ) = +65.0° 
and £i7r( 3 Pi) = —34.6°. If one considers all those one-loop box diagrams which include in the 
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upper part all tree graphs evaluated in section 4.2 and in the lower part only the one-pion 
exchange one gets for the imaginary part of £>, 

ImB = B (thy) • M 3 Pi) = 2.74 fm 4 • (-0.604) = -1.65 fm 4 , (B.5) 

a value which is in 10% agreement with the empirical Imi3 = —1.5 fm 4 . Of course, this merely 
reflects the fact that the empirical 3 Pi phase shift at is not far from one obtained in lir- 
exchange approximation. Note, however, that for the 3 Po phase shift at the empirical 
value and the l7r-exchange approximation differ by a large factor of —10. 



C Effective field theory approach to final— state inter- 
action 

In this appendix we want to give an elementary derivation of the final-state interaction cor- 
rection factor Eq.(6), F p (W) = [1 + a 2 P 2 ]^, P 2 = W 2 /A - M 2 , in scattering length ap- 
proximation (i.e. for r p = 0). Close to threshold all final state three-momenta are small 
and therefore one can approximate both the meson production process NN — > NNtc and 
elastic iViV — > NN scattering by momentum independent contact vertices proportional to A 
and the scattering length a p , respectively. Consider first low energy iViV-scattering in this 
approximation. The bubble diagrams with 0,1,2,. . . rescatterings can be easily summed up in 
the form of a geometric series, 

a ^^\ (2 7 r)3^ + ^-P + - = flp + t ^ + - = T^:' (CU) 

using dimensional regularization to evaluate the (vanishing) real part of the loop integral. 
Obviously, the sum of these infinitely many loop diagrams is just the unitarized scattering 
length approximation which leads to, 

tan 5 (W) = a p P* . (C.2) 

Next, consider in the same approximation meson production followed by an arbitrary number 
of iViV-rescatterings in the final-state. Again, these loop diagrams can be summed up to, 

A (C3) 



and taking the absolute square, 



A 



1 (XpPjf, 



1,4,2 (C.4) 



1 + ajP? ' 



one encounters the final-state interaction correction factor F P (W) = [l + a p (W 2 /4 — M 2 )] -1 in 
scattering length approximation. Since, the scattering length is much bigger than the effective 
range parameter for iViV-scattering (a p » r p ) one has already derived the dominant effect 
due to final-state interaction. Of course, in order to be more accurate one should eventually 
go beyond momentum independent contact vertices. The main point we want to make here is 
that the final-state interaction correction factor F P {W) (for r p = 0) has a sound foundation 
in effective field theory. 
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Figures 




Figure 1: The process pp — > ppir in the center-of-mass system. 
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Figure 2: Fit to the total cross section for pp — » pjwr as described in 
the text (solid line). The data are from (boxes) and @ (crosses). The 
dashed line is explained in sec. 3.4. 
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Figure 3: Feynman graphs for neutral pion production, (a) and (b) 
are the direct terms, with meson exchange M. = 7r°, 77, u, p°. (c) is the 
rescattering graph. The heavy dot denotes the second order isoscalar chi- 
ral 7r7riViV-vertex. Graphs where the pion (dashed line) is emitted from 
the other proton (solid) line and graphs with crossed outgoing proton 
lines are not shown. 
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Figure 4: Class of one loop graphs involving the nn interaction. For 
further notation, see Fig.3. 
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Figure 5: Further loop diagrams proportional to <jvjv/ ft- ( a ) an( i (b) 
represent 7r + 7r _ exchange between protons, (c) and (d) are the n + rescat- 
tering diagrams. For further notation, see Fig.3. 
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Figure 6: Loop diagrams with ir° emission via charge exchange. For 
further notation, see Fig.3. 
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1875.0 1877.5 1880.0 1882.5 1885.0 1887.5 1890.0 

W [MeV] 

Figure 7: Pion production cross section as a function of the pn invari- 
ant mass W . The three curves corresponds to T lab = 299.3, 306.3 and 
314.1 MeV, in ascending order. 




Figure 8: The eta-production cross section a tot (pp — > ppr/) as a function 
of Tiab. The data are taken from ref.0. 
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Figure 9: Diagrams that give rise to the non-vanishing imaginary parts 
Im.4 and ImB in one-pion loop approximation. The projection on the 
on-shell pp intermediate state is symbolized by the double line cutting 
the diagrams. For further notation, see Fig. 3. 
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